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Abstract — Cascading failure of a power transmission system 
are initiated by an exogenous event that disable a set of elements 
(e.g., lines) followed by a sequence of interrelated failures (or 
more precisely, trips) of overloaded elements caused by the 
combination of physics of power flows in the changed system 
topology, and controls. Should this sequence accelerate it can 
lead to a large system failure with significant loss of load. 
In previous work we have analyzed deterministic algorithms 
that in an online fashion (i.e., responding to observed data) 
selectively shed load so as to minimize the amount of lost load 
at termination of the cascade. In this work we present a rigorous 
methodology for incorporating noise and model errors, based on 
the Sample Average Approximation methodology for stochastic 
optimization. 

I. Introduction 

We present a rigorous methodology for computing ro¬ 
bust algorithms for control of cascading failures of power 
transmission systems. We focus on the linearized, or DC, 
approximation for power flows, and on line tripping outages. 
We note, however, that the same underlying methodology 
will apply for other models of power flows and other types 
of equipment outages. This work extends the approach in 
[3]; related work is consided in [7], [4], [8], and references 
therein. 

In the DC approximation to power flows (see [2] for 
background) a transmission system is modeled by given a 
directed graph G with n buses and m lines. In addition, for 
each line e we are given its reactance x e and its limit u e (we 
may also refer to a line e in the form e = pq so as to indicate 
its “from” bus p and ‘its ‘to” buse q. Additionally, we are 
given a supply-demand vector j3 £ lZ n with the following 
interpretation. For a bus i, if 3, > 0 then * is a generator (a 
source node) while if /3j < 0 then * is a load (a demand node) 
and in that case — pi is the demand di at i. The condition 
Bi = 0 is assumed to hold. We denote by Q the set of 
generators and by V the set of demand nodes. The linearized 
power flow problem specifies a variable f pq associated with 
each line pq (active power flow) and a variable cf) p (phase 
angle) associated with each bus p. The DC approximation is 
given by the system system of equations: 

Nf = /3, N T (/) — Xf = 0, (1) 
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where N denotes the node-arc incidence matrix of G [1] and 
X = diag{a,’ij}. 

Remark 1.1: It can easily be shown that system 0 is 
feasible if and only if YlieK A; = 0 for each component 
(“island”) K of G, and in that case the solution is unique in 
the / variables. 

In this statement the line limits u pq play no role. It may 
be the case that in the (unique) solution f to 0 we have 
\fpq\ > u pq . In that case line pq is at risk, and, unless some 
control action is taken, will eventually trip. The “eventually” 
is an imprecise statement. Certainly, in the context of cascade 
modeling, the u pq should be the “emergency” line limits. 
Nevertheless, care should be taken to model the non-precise 
nature of line tripping. We will comment on this point, again, 
later. 

We can now outline the model for discrete-time cascading 
failures from [3]. Similar models have been proposed in the 
literature, see 


Template 1: GENERIC CASCADE TEMPLATE 
Input: a power grid with graph G (post-initiating 
event). Set G 1 = G. 

For t = 1,2,... Do 

(comment: time-step t of the cascade) 

1. Set /* = vector of power flows in G*. 

2 . Set O t = set of lines of G 1 that become outaged 

in time-step t. 

3 . Set G t+1 = G t — O t . Adjust loads and 

generation in G f . 

In this template, each “time-step” models a time increment 
of length At > 0. In Step 3, when a load imbalance results 
in some island, we assume that the larger of generation 
or load are adjusted, downwards, so as to attain balance. 
The model does not attempt to explain the cascade within 
resolution finer than At. It should be noted that cascades 
can be extremely rapid - perhaps lasting seconds. However, 
several cascading failures of notoriety have been significantly 
slower, such as the 2004 U.S. Northeast event [5], For the 
purposes of this paper, the reader should assume that At is of 
the order of minutes. We remark that cascades are extremely 
noisy events, characterized at a fine time scale by a myriad of 
events, some of which are extremely challenging to model, 
such as physical contacts of sagging lines with vegetation, 
and human action (and errors - which have been found 
to take place in many cascades). Any model of cascading 
failures which attempts to model time at a very fine scale 
(perhaps, continuous time) would need incorporate corre- 



spondingly fine models for such events. Additionally, there 
would be a need to resolve “race-conditions,” which take 
place when many lines can trip within a very short amount 
of time with potentially very different cascade outcomes 
(observed by other authors). In this paper, for simplicity, we 
will therefore assume a relatively large At. We note that 
one could furthermore argue that a control algorithm that is 
subject to noise and incomplete information would benefit 
from relying on a relatively large At. 

In [3] we considered a general form of control template, 
which assumes that the initial exogenous event that sets- 
off the cascade has been observed, and that G is the post¬ 
event network. The control entails load shedding; when load 
is shed in some connected subnetwork, generation must be 
correspondingly adjusted in that subnetwork. 

Framework 1 Cascade Control_ 

Input: a power grid with graph G. Set G 1 * 11 = G. 

1 : Compute control algorithm 

2 : For f = l,2,...,T—1 > controlled time-step t. of 

the cascade 

3: Obtain grid measurements 

4 : Apply control 

5: Set /W = vector of resulting power flows in 

6 : Set 0 W = set of lines of that become outaged 

in time-step t 

7: Set G^ t+1) = GW - 

8 : Adjust loads and generation in G® 

9: Termination (time-step T). If G rr> has line overloads, 
proportionally shed demand until all line overloads are 

eliminated. Set T "l = min 11, maxj { ^ - } j. If 
ipt 7 " 1 > 1, then any bus v resets its demand to /ip^ T \ 


We stress that G is the post-initiation state of the grid, 
i.e. the elements disabled in the initial exogenous event are 
not present. The termination condition ensures that indeed 
the cascade is stopped at the end of the planning horizon 
[3]. In general. Step 3 consists of any step where a network 
controller obtains current data. In Step 4 this data is used to 
apply the control computed in Step 1 using the measurements 
as inputs. A completely deterministic version of Framework 
[T] is considered in [3]. This algorithm is characterized by two 
features: 

(s.i) At time t, if I\ is an island in G (4 \ then in Step 4 all 
loads in K are scaled by a factor 0 < \ f K < 1 (as are 
all generator outputs in K). 

(s.ii) In Step 6, a line e is tripped if u e < lsi 4 ) l- 
The objective is to compute the A 4 so that at termination 
the total load still being delivered is maximized. The main 
contribution in [3] is an algorithm that solves this problem 
in polynomial time, for each fixed value of T. This fact may 
appear surprising given the nature of rule (s.i). There are 
an exponential number of islands K that could be realized 
at time-step t - how could the algorithm run in polynomial 
time? 


The answer is that the A* parameters need only be com¬ 
puted for an optimal realization of the controlled cascade. 
Having chosen the vectors A 1 , A 2 ,..., A 4-1 then a unique 
state will be observed in time-step t. This is a consequence 
of the fact that rule (s.ii) is purely deterministic. Hence if we 
somehow know that A 1 , A 2 ,..., A 4-1 have been optimally 
chosen, then all that is needed is for A 4 to be optimally 
chosen, as well. Continuing inductively we will obtain an 
optimal control across all time-steps. To put it in a pedestrian 
manner, the controlled cascade plays out like a script, with 
the events that transpire at each time-step t known, precisely, 
in advance. 

Thus, the key in the analysis of the algorithm in [3] 
is that, indeed we can choose A 4 optimally (provided that 
A 1 , A 2 ,..., A 4-1 have previously been computed optimally). 
To prove this point [3] relies on a variant of dynamic 
programming. Below we will consider an updated form of 
this algorithm. 

II. Modeling stochastics 

Even though the algorithm in [3] is provably optimal, it 
is readily apparent how that algorithm falls short, and in 
particular may not prove robust. This concerns rule (s.ii) - 
should this assumption prove inaccurate it is quite likely that 
the set of lines that trip in Step 6 will be different than 
anticipated in the “script” mentioned above. Not only that, 
but the set of islands actually observed at time-step t may 
be different from those expected by the “script”, and thus 
the computed control does not even make sense (i.e. we 
have the wrong parameters A^-). In this section we consider 
algorithms that not only bypass these shortcomings, but also 
attain a form of algorithmic robustness that can be precisely 
stated. 

(L.a) We will incorporate stochastics into the line-tripping 
rule. More precisely, at time t we compute, for each 
line i, the quantity ff = (1 + e\)\f^ |, where where 
e 4 is a random variable from a known distribution. 
Additionally, we assume that |e 4 | is bounded by 1 > 
b > 0. Here, it is assumed that b « 1. The stochastic 
line-tripping rule is that line i is tripped if / 4 > it,;. In 
what follows, we will / 4 is the noisy flow in line i at 
time-step t. 

(L.b) Under the stochastic line tripping rule, control rule 
(s.i) does not make sense because we do not know, 
in advance, the set of islands that will be realized at 
time-step t. In fact, potentially, any island could be 
realized. Instead, the control will compute a single value 
A 4 which will be used to scale all loads at time-step t. 
We can now state the optimization problem of interest. 

Definition 2.1: Optimal Robust Control Problem. Com¬ 
pute scaling values A 1 ,..., A T_1 such that subject to rules 
(L.a), (L.b) we maximize the expected yield , where by 
“yield” we mean the load that is served at termination. 

Remarks. 

• In (L.a) we do not make any assumptions as to the 
source or nature of the stochastics. In fact, we gener- 
ically model “noise” in this fashion, so as to be able 








to capture any form of uncertainty that could hamper 
a control algorithm, so long as the magnitude of the 
errors is not overly large. As we will see in our 
experiments, simply allowing stochastic tripping gives 
rise to a large variety of cascading outcomes. A control 
that maximizes expected load will thus be robust with 
respect to many alternative histories that the system 
could follow. However, there is a specific setting in 
which rule makes sense - we can use (L.a) to measure 
errors in line measurements. 

• The termination rule in Framework[[]needs to be revised 
when rule (L.a) is applied so as to make certain that the 
final load shedding does terminate the cascade. We do 

f /- T> 1 

so by redefining = min < 1, maxj{ (1 ~b)uj } f 

Since b is small, this amounts to a small correction, at 
termination. 

III. Solving the optimal robust control problem 

UNDER A GIVEN NOISE VECTOR £ 

Some of the algorithmic steps presented here echo some 
steps in [3]. However, because of (L.a) and (L.b), the under¬ 
lying mathematical nature of the problem is fundamentally 
different, as are the actual proofs. First we introduce some 
notation. 

• We define the function as the maximum 

expected yield given that in the first time-step the noisy 
flows are /W = zf. 

• We denote G 1 the power grid at the beginning of time- 
step t, a random variable under (L.a). 

• Denote by e the m x T matrix of f'f 1 values for all i 
and 1 < t < m. 

• For z > 0 real, let @^(z|/, e) be the deterministic 
maximum yield obtained if at the start of the first time- 
step the flows in the grid are equal zf, and all the e is a 

given, fixed, vector of values e'p (rather than random). 

(T) 

Note that 0) ; ’ is deterministic since all realizations of 
the random variables ep are provided. 

Clearly the following holds, 

4 T) w/~) = E e 

where /' = 

Lemma 1: e) is a nondecreasing piecewise- 

linear function of z with two pieces, the second one of which 
has zero slope. 

Proof: Note that since T = 1, only termination step 
in Framework [l] will be executed. Denoting by D the sum 
of demands implied by /' we have as per our cascade 
termination criterion that the final total demand at the end 
of T = 1 time-step will equal 

zD , if z < 1 /t/> (1) , and (3) 

Wb D otherwise - (4) 


Using similar ideas to those shown in [3], we now turn to 
the general case with T > 1 and we will show the following 
theorem. 

Theorem 1: Qp\z\f, e) has at most ( m br +1 y break¬ 
points. 

To prove Theorem [T] we will need the following defini¬ 
tions. 

Definition 3.1: A critical point is a real 7 > 0, such that 
for some line j, 7 fj = Uj, i.e. y\fj\ = Uj/( 1 + e‘). 
Remark: for completeness we should use, in this definition, 
the superindex t, which we have skipped for brevity. Recall 
that we assume Uj > 0 for all j; thus let 0 < 71 < 72 < 
... < 7 P be the set of all distinct critical points. Here 0 < 
p < m. Write 70 = 0 and y p +i = + 00 . 

Definition 3.2: For 1 < i < p let F W = {j : 7 h fj = 

Uj}. 

Assume that the initial flow is zf with 1 > z > 0 and 
let 0 < A ' 1 > < 1 be the optimal multiplier used to scale 
demands in time-step 1. Write 

q(z) = argmaxj/i : 7 h < z}. (5) 

Thus, z < 7 q +i, and so A^z < 7 g +i. We stress that these 
relationships remain valid in the boundary cases q = 0 and 

q=p. 

Notation 3.1: Let the index i be such that A^z £ 

( 7 i-i, 7 i]- 

In time-step 1 of Framework [T] at line [4] we will scale 
all demands by A^ 1 7 We assume the framework is given 
a connected graph, namely, G^ l) is connected. Therefore, 
in line [ 5 ] we will also scale all supplies by A' 1 7 Thus, for 
any h < i — 1, and any line j £ F^ h \ we have that after 
Step [ 5 ] the absolute value of fj is X^zfj > 7 hfj = Uj, 
and consequently line j becomes outaged in time-step 1. On 
the other hand, for any line j ^ L>h<i-iF^ h ], the absolute 
value of the flow on j immediately after Step pi is \ < ] Gfj < 
jifj < Uj, and so line j does not become outaged in step 1. 
In summary, the set of outaged lines is in other 

words, we obtain the same network G^ = G^ 1 pU t h p 1 F^ 
for every z with A^z £ ( 7 »-i, 7 »]. 

Notation 3.2: For an index j, write lC(j) = set of compo¬ 
nents of \ U J h=1 F( h \ 

Let H £ K,[i — 1) and denote the initial supply-demand 
vector by 77 Prior to line [ 8 ] in step 1, the supply-demand 
vector for H is precisely the restriction of \ <r> z3 to the 
buses of H, and when we adjust supplies and demands in 
line [ 8 ] we will proceed as follows 

• if E se x>rw(- A(1) ^) > 'ZsegnH( x( 1 )z Ps) then for 
each demand bus s £ 2? fl H we will reset its demand 

to —rX^zfis, where 

_ E s ggng( A ^ z ^) _ _ Esggo ipJppf 
Yhs&DC\H^~ X ^ z Ps) EsGX>n h(P s ) 

and we will leave all supplies in H unchanged. 

. likewise, if E se cn h(~ zX< " 1) Pb < 

E s ggnt?( A(1 ^^ s ) ti 1611 tf 16 supply at each bus 




( 2 ) 










s £ Q D H will be reset to rX 1 ^zB s , where 

r= _T^nM Wsl 

X^segnH 

but we will leave all demands in H unchanged. 

Note that in either case, in time-step 2 component H will 
have a supply-demand vector of the form X^z/3 H , where 
is a supply-demand vector which is independent of z. 
The supply-demand vector corresponds with flows f 11 
on the lines in H, which are therefore also independent of z. 
The flows /^ on G ( 2 ) 

are U ff 6 G( 2 )/ ff , and the final total 

demand will be 

eg-^AWzi/W) ( 6 ) 

For a given value of i, since z < 1, also i < q holds. As 
noted above, by definition ([5]) of q we have that 7 , < z q < 2 . 
Thus, the expression in (| 6 p is maximized when A f l * = 7i , 
and we obtain final (T-step) demand equal to 

A = u ff 6 jc(i-i) f H ), (7) 

In summary, 

®Gm(. z \f) = max A (8) 

We are now ready to prove Theorem [I] 

Proof: In Lemma [T] we showed that for T = 1 the 
function has two breakpoints. As can be seen from (| 8 j, the 
number of breakpoints is at most q(z). It is important to note 
that (a) q(z) is at most the number of lines in the graph and 
(b) the number of non-tripped lines decreases by at least 1 in 
every time-step (otherwise, the cascade stopped). Therefore, 
the number of breakpoints for T steps is m(m — 1)... (to — 
T +1)=(F^t)T- ■ 

IV. Solving the optimal robust control problem 
using the Sample Average Approximation Method 

The key methodology that we will reply in order to solve 
the problem of interest is the Sample Average Approximation 
(SAA) method [ 6 ]. The method begins by taking mTN i.i.d 
samples (each under the assumed distribution) of e values. 
We view these values as arranged into N ensembles, or 
realizations. Here, for 1 < k < N we let e k denote the 
k th realization, consisting of values e^’ k , for 1 < t < T 
and all lines i. 

Definition 4.1: Given any candidate vector A = 
(A^ 1 ), A^ 2 ),..., A( T-1 )) for solving the robust optimal 
control problem, denote by 0g(z|/, A, k) to be the yield 
obtained at termination if at the start of the first time-step 
the flows in the grid are equal zf. if we use A as the control 
vector at each time-step, and is used for each t and i 
in the line-tripping rule (L.a). 

Thus, ©g(1|/ 1 i A, k) indicates the behavior of the con¬ 
trol vector A under the fixed noise vector e k . This motivates 
the following definition. 


Definition 4.2: (Sample average robust control prob¬ 
lem.) Given the N realizations e 1 ,... ,e N , compute a control 
vector A so as to maximize 

1 N 

^Eesai/N.A,*). (9) 

v fc=1 

This definition is appealing in that, clearly, if N is large a 
control that maximimized ([9]) should clearly perform well in 
our stochastic setting - since it maximizes the average out¬ 
come over many possible noise scenarios. This observation 
can in fact be rigorously established by using the Central 
Limit Theorem (and observing that all random variables 
under consideration are bounded). See e.g. equation (2.23) 
from [ 6 ]. If we want to solve the Robust Optimal Control 
problem within additive error <5 > 0 with probability at 
least 1 — a one can establish an upper bound of the form 
N = log(l/a), where o is the standard deviation of 

the distribution of the e values [j] We can establish an analogue 
of Lemma|T] Let rfc{z\f) = ± J2k=i 0 g(“I/> A > k )- This 
is the maximum average yield at termination under all the 
realizations if we start time-step 1 with flows zf. 

Lemma 2: fjc{z\f) is nondecreasing, piecewise-linear, 
with at most 

N ——- 

(to - T)\ 

breakpoints. 

Proof: By theorem [TJ Q^\z\f, e) is piecewise-linear 
with at most breakpoints. Since f]Q(z\f) is an 

average taken over N instances of Q < o'\z\f,e), the lemma 
follows. ■ 

Using Lemma [2j it follows that we can efficiently deter¬ 
mine the control in order to maximize the expected total 
demand in the face of measurement errors. 

V. Simulation Study 

In this section we present simulation results using a mod¬ 
ified form of the IEEE 118-bus system, with line limits set 
approximately at 20% above flow values^] We first describe 
our implementation and specific noise model (see (L.a) 
above). We then present and discuss the computed i)g(z|/), 
obtained via the SAA, and compare the performance of our 
robust control algorithm to its non-robust version from [3]. 

A. Implementation and Noise Model 

As discussed above an approximately optimal load shed¬ 
ding control can be found by computing the function 
flj.{z\f), and for this purpose we will rely on Lemma [I] 
Thus, at each every time-step t, we solve the DC equations 
to obtain the flows and compute the critical points 71 ,..., / y p . 
An important point is that these are the critical points arising 
from all realizations (recall Definition [3TT|. Next, for each i 
with 1 < i < p we proceed as follows. 

Note that any two choices for A 1 in the open interval 
( 7 i_i, 7 i) will result in exactly the same set of line trips 

1 Some technical points skipped for brevity. 

- We will present experiments with larger examples at the conference 







(a) (b) 


Fig. 1. The yield T = 3 in the first time step for: (a) 5 realizations, and 
(b) 120 realizations. 



(a) (b) 

Fig. 2. The yield when considering 120 realizations for: (a) the second 
time step and for T = 3, given that an optimal load shedding scaling is 
done in the first time step, and (b) the first time step and for T = 2. 


under all realizations. In fact, any two choices for A 1 in 
( 7 i_i, 7 i) will produce the same graph G{i ) (which we can 
easily compute) but with supply-demand vectors that differ 
by a constant factor - and thus, in a common flow vector 
(up to scale) which we denote by f(i). It follows that if we 
(recursively) compute the function Vop{z \/(*)) (which is a 
function of the real z) we will obtain flQ{z\f) in the interval 
(7*-i> Til- 

Next we discuss the specific noise model used to imple¬ 
ment rule (L.a). We will first describe a specific noise model. 
Then we will prove it is of the form (L.a). A line j will trip 


with probability 1, if fj > Uj (10a) 

with probability 1/2, if Uj > \fj\ > 0.95 Uj (10b) 
with probability 0, otherwise. (10c) 


Lemma 3: Rule (lOi is an example of (L.a). 


Proof. Define the random variable e 
follows: 


(t) 


for every j and t as 


jt) _ 0 05 ^ with probability 0.5, and (11) 

J 0.95 

e'p = 0, with probability 0.5. (12) 

Recall that ftp = (1 + e'p)fj and consider the case 
where fj > 0.95it,,-. With probability 0.5, cp = 0 and 
therefore fj 0) < u ;) and line j will not be tripped. But, with 
probability 0.5 we get r { p = and /jW = (1+ [/§§)/■ > 

(1 + 0.95m.? = 11 j an d li ne 3 will trip- ■ 


B. Results 

We tested values of T = 2,3,4, 5. The initial cascade is 
caused by tripping line (4,5), since tripping this line causes 
a relatively major cascade in the 118-bus system. 

Fig. □ shows the yield vs the scaling in the first time 
step for a value of T = 3. In Fig. OT) the number of 
realizations used for the Sample Average Approximation 
(SAA) is 5 while in Fig. |TJb) this number is increased to 120. 
As expected, increasing the number of realizations affects the 
expected yield in the different linear segments. In Fig. |T| we 
see that despite the fact that the graph differ, in this specific 
case the maximum yield remains similar. 



Fig. 3. The yield when considering 120 realizations and T = 5 for: (a) 
the first time step and for T = 2, and (b) the second time step, given that 
optimal load shedding is performed in the first time-step. 


While Fig. [I] shows the load shedding required in the 
first timestep (shedding of « 0.72 maximize the yield), the 
figure does not reveal what will be the load shedding at the 
second time step. Fig. |2|a) shows the yield as a function of 
load shedding in the second times tep assuming that in the 
first time step an optimal choice for A 1 is used. As can be 
observed from Fig. [2ja), the optimal choice for A 1 is ss 0.85. 
Fig.[2jb) considers the case where 7' = 2, it can be seen that 
that load shedding is done at only one time-step. While the 
maximum yield is slightly lower, overall it is close to the 
maximum yield (« 62%) for the case where T = 3. We will 
now show that the difference in the maximum resulting from 
increasing T can be significant. 

In Fig. [3] we set T = 5. Fig. [3ja) shows the yield as a 
function of load shedding in the first time steps and can be 
used to find the optimal load shedding scaling for the first 
time step. Fig. [3jb) shows the yield as a function of load 
shedding in the second time step given that the optimal load 
shedding scaling was done in the first time step. We can see 
that in this case making a small amount of load shedding in 
the first two time steps results with a very good final yield 
(~ 85%) compared to the case of T = 3. We remark that in 
the rest of the time steps the optimal control is to carry out 
very minimal load shedding. 

A natural question is what would the non-robust control 
algorithm from [3] would obtain in the considered model. 
To answer this question, we compute the control using the 
algorithm from [3]. As expected, in the non-robust model, 
the algorithm from [3] is the best, since there is perfect 









T 

2 

3 

4 

5 

Non-robust solution 
and non-robust model 

65.46% 

65.46% 

74.44% 

86.84% 

Non-robust solution 
and robust model 

31.92% 

30.46% 

47.75% 

23.07% 

robust solution 
and robust model 

62.19% 

62.19% 

70.73% 

78.36% 


TABLE I 

Yields obtained from performing the solution obtained by 

THE ROBUST AND NON-ROBUST ALGORITHMS UNDER THE ROBUST AND 
NON-ROBUST MODEL. 


T 

2 

5 

5 realizations 

[62.39%, 63.19%] 

[83.02%, 84.87%] 

60 realizations 

[62.7%, 63.08%] 

[83.84%, 84.62%] 

120 realizations 

[62.79%, 63.04%] 

[84.05%, 84.64%] 


TABLE II 

A 95% CONFIDENCE INTERVAL FOR OPTIMAL EXPECTED YIELD FOR 
DIFFERENT NUMBER OF REALIZATIONS. 


knowledge of the flows and the algorithm from [3] is 
optimal. However, when measurement errors are present, the 
algorithm from [3] performs much worse than the robust 
algorithm presented here. Table [I] shows the yield obtained 
by both algorithm. Additionally, Table [I] further demonstrates 
that increasing T improves the yield. We remark that the non- 
robust model refers to the model from [3] and the robust 
model refers to the model considered in this paper. 

We also compute a 95% confidence interval on the yield 
obtained when performing the optimal load shedding. The 
confidence interval is computed for different number of re¬ 
alizations and T = 5 and it is shown in Table HT1 The obtained 
confidence interval are overall very good. As expected, we 
can see improvement in the confidence interval as we take 
more realizations. 
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